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Abstract. Disease awareness in epidemiology can be modelled with adaptive 
contact networks, where the interplay of disease dynamics and network alter- 
ation often adds new phases to the standard models ((T|, (2)) and, in stochastic 
simulations, lets network topology settle down to a steady state that can be static 
(in the frozen phase) or dynamic (in the endemic phase). We show for the SIS 
model that, in the endemic phase, this steady state does not depend on the ini- 
tial network topology, only on the disease and rewiring parameters and on the 
link density of the network, which is conserved. We give an analytic description 
of the structure of this co-evolving network of infection through its steady-state 
degree distribution. 

1 Elaborating on the Pairwise SIS Model with Rewiring 

Gross et aL (fT)) incorporated disease awareness into conventional SIS dynamics by al- 
lowing susceptibles to evade infection through retracting links from infected neighbours 
and rewiring them to randomly selected susceptibles. Using the moment closure from 
O to obtain a pair approximation model (PA) of the process, the time evolution of the 
fraction of infected i is coupled to that of the densities Zxy of links connecting nodes 
of type X and Y ({X,Y} G {S= susceptible, /= infected}). Apart from the con- 
served link density I ~ Iss + hi + hi, the model can be parameterised by p, w S [0, 1], 
with rewiring, infection and recovery rates given by oj, {1 — uj) p and (1 — w) (1 — p), 
respectively. 

The full phase diagram is presented in Fig.[Tl for any ^ > 0.5 the minimum rewiring 
rate for bistability to arise is a; = 0.5. The transcritical bifurcation, serving as the epi- 
demic threshold and approximated in 1 1], occurs at exactly oj = 1 — 1/ {p[2l + 1]). 

We have checked the steady state values for densities and pairs predicted by the PA 
against simulations on networks with N = 10^ nodes, taking as initial conditions the 
predicted steady-state value of i, with infected and susceptibles initially randomly dis- 
tributed on a Poisson network of average degree (k) =21 = 20. In the simple endemic 
phase, prevalence and link densities of numerical simulations and model predictions 
agree reasonably well. At w = p/ (1 + p), i.e. when the infection rate is equal to the 
rewiring rate, the pair approximation model correctly predicts i and the mean degrees 
of both susceptibles and infected. 
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Fig. 1: Phase diagram of pairwise equations in [I] , / — 10. Compiled widi (3- 

2 Observing A Steady- State Topology 

The steady-state degree distributions (DD) of the subnetworks of infected and suscepti- 
bles - as well as that of the complete network - are shown in Figs.|2]and[3] Independently 
from topological starting conditions in the simple endemic phase, the final topology (as 
described by several standard measures) settles down to a characteristic state (Fig.|2l). 

Let (fcs) = J2k k^k/s and (fci) — J2k ^^k/i be the average degree of susceptible 
and infected, where is the fraction of the total number of nodes which are in state 
a and have degree k, a G {s, i}. Hence the average degree of the network is (fc) = 
s{ks) + i{ki). We have also, by definition, s(fcs) — 2/ss + 'si and i(fci) = 2lu + /si, 
so that the average degrees of the two subnetworks in the steady state are determined 
by the equilibrium values of the PA in the region where the moment closure from li3] 
holds. 

Rewiring induces infection-status clustering (|1|); in the PA, the difference of the 
mean degrees of susceptible and infected nodes is independent of the link density and 
given by 

(1 - w)p 

Thus only for uj > pj (l + p) does rewiring yield a bigger mean degree in susceptibles 
than in infected (Fig.O. Again uj and p given by a; = p/ (1 + p) enjoy a special status, 
in that then the mean degrees of the S- and /-ensemble are the same, as are their degree 
distributions and every other topology -related statistics we measured. 




Fig. 2: Initial topologies settling down 
to final DD. Erdos-Renyi (ER), Watts- 
Strogatz with rewiiing probability /3 — 
0.1 (WS), Barabasi-Albert (BA). N = 
10^, (k) = 20,LO = 0.9, p = 0.7, 
averaged over 10"^ samples (between 
4000 < t < 5000 for final DD). 




(a)a; = 0.1 (h)ui = 0.9 

Fig.3: Infection-status clustering in the steady state in stochastic simulations (MC). Fig. [3a] A{k)MC = 

-0.50, Zi(fc)pA = -0.47. Fig.|3b] z\{fe)MC = 8.01, z\(fe)pA = 6.35. Initial Poisson DD, AT = 10^ , (k) = 
20, p = 0.7, averaged over 10* samples between 1000 < t < 2000. 



3 The Node Cycle 

Depending on a node's status, it will either gain in degree (when susceptible) or lose 
links (when infected). Characterising such a typical node cycle gives important single- 
node characteristics describing its degree during and lifetime of each stage, as well as a 
means to determine the overall infection-class degree distributions. 

The time evolution of the number of susceptible and infected neighbours Is (t) and 
li{t) of a typical node is 

lsi{t) =cie-^i* + A2i + C2 

ln{t)=C3e-^'' + X3t + C4 (1) 
in its susceptible stage (phase I) and 

/S2(i) =C5e-^**+C6e-^=* 

luit) = C7e-^^' + cse-^'' (2) 



when it is infected (phase II). The constants A,; e M>o and Cj e R are comprised of 
w, p, the steady-state values of /gi and s (computable through the PA), as well as the 
four initial conditions to be determined by a cyclic closure of ^ and (|2]l. 

A crucial ingredient is the mean-field approximation for the average lifetime rs 
of a susceptible; ts = s/ ([1 — w] p^si) turns out to be a very good one throughout 
parameter space. This is due to the fact that over a susceptible's lifetime, there is no 
noticable degree-, time- or any other correlation in the 5-nodes rewired to it, particularly 
not in the number of their /-neighbours. Hence the node-cycle model also accurately 
describes regimes with very high lj. There, the near-linear degree growth and the trough 
of li{t) around rs in the first phase are well-captured by our equations. So is the rapid 
rewiring away of /s -links in the initial stage of the second phase, followed by slow 
degree loss through recovery and subsequent rewiring (Fig.|4aland|4b]l. 




(a) Model (b) Simulations 

Fig. 4: Degree evolution of a single node, w = 0.9, p = 0.7, / = 10. Fig.|4a] isi (0) = 0.43. ;ii(0) = 15.6, is2(0) = 
isi(Ts),ii2(0) = iii(Ts).Fig.|4b]A' = 10''. 



4 Deriving the Degree Distribution 

Let Zsi^ be the fraction of the S'/-link density with /-nodes of degree k. Denote also 
by /sfci the number of 5/-links whose 5-node is of degree fc, divided by the number 
of nodes N . In order to describe the infection-class subnetworks, we must, instead of 
considering the pairwise equations in fTj, deal with the dynamics of each degree class of 
infected and susceptibles. For the processes and rates we are considering, this is given 
by 

- (1 - w)(l - p)ik + (1 - w)pls^i + w (Zsi,+i - ku) 
{1 - w){l - p)ik - {I - w)pls^i+ w— {sk-i - Sk) ■ (3) 



dik 
dt 
dsk 
~dt 



4.1 Inside the PA regime 



Assuming that the end nodes of ^si-links uniformly sample the two subpopulations, (O 
may be closed by taking 

7 7 ^■Sfe A 7 7 

From ([3]) with i, s, (/cg) and (fci) in (|4]i given by their values in the steady state of the PA 
in [1 1, we obtain recurrence relations for the steady-state values of and i^. Together 
with the normalisation conditions, these equations can be solved for the two subnet- 
work's stationary DD. In the region where the PA performs well, this is a reasonable 
approximation (results for a sample point in this region are shown in Fig. |5a] Also in 
this region, state-degree correlations are negligible, and the DDs of the two subnetworks 
are practically identical (as observed in simulations before). 



4.2 Outside the PA regime 

The fully analytic, PA-based procedure described above breaks down whenever the PA 
does, and state-degree correlations then become important. Since the PA no longer holds 
for the (fcg.i), these cannot be treated as parameters anymore and must be expressed in 
terms of the variables Sk, ik in ©. Because of the state-degree correlations, a different 
closure assumption than (|4]i must be used for Instead of the second equation in (|4]l, 
we shall take 

where a G [0, 1] is a phenomenological coefficient obtained through least squares fit- 
ting. 

The fact that a must be nonzero can be understood as follows: Links of type SI 
are very short-lived with respect to the average infectious period. This means that the 
/-nodes of 5'/-links do not uniformly sample the infected population, and that the sam- 
pling is biased towards the 'younger' infected nodes, born from susceptible nodes. This 
sampling bias can then be modelled as a weighted average of the Sk- and u-distribution. 
This argument holds throughout, but only when state-degree correlations are important, 
i.e. for w > yo/ (1 + p), is this correction important. A similar argument applies to the 
ZSfci-distribution, but the sampling bias is far less pronounced in this case, because the 
average susceptible period is smaller than the average infectious period in this region. 

With these two modifications with respect to the preceding section, the equilibrium 
of ^ can also be found (although no longer through recurrence relations) and compared 
with the stationary degree distributions obtained from simulations (Fig. l5b] i. The best 
choice of the coefficient a gives a good quantitative agreement for the whole range of 
the DDs; whereas the predictions of the naive PA-based approach used before are way 
off the Monte Carlo results. 




(a) uj = 0.5 (b) cj = 0.9 

Fig. 5: Analytic stationary DD of infected and susceptible subnetworks vs. MC results. Simulations with initial Poisson DD, 
Af = 10*, (fc) = 20, p = 0.7, averaged over 10"* samples between 1000 < t < 2000. 
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